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ABSTRACT 

The GeV light curve of a pulsar is an important probe to detect acceleration regions in its magnetosphere. 
Motivated by the recent reports on the observations of pulsars by Fermi Large Area Telescope (LAT), we 
restudy the two-pole caustic model and revise it to investigate the properties of the light curves in the GeV band. 
In the revised model, although acceleration gaps can extend from the star surface to the light cylinder along near 
the last open field lines, the extension of the gaps along the azimuthal direction is limited because of photon- 
photon pair production process. In such gaps, high-energy photons are emitted uniformly and tangentially to 
the field lines but cannot be efficiently produced along these field lines where the distances to the null charge 
surface are larger than ^ 0.9 times of the distance of the light cylinder, and the effective azimuth extension of 
the gaps is about 230°. The model is applied to the four pulsars Vela, PSR J1028-5819, PSR J0205H-6449, and 
PSR J2021H-3651 whose light curves obtained with Fermi have been recently released. The model is successful 
in reproducing the general feature of the light curves for the four pulsars, and the radial distances of the radio 
pulse for the four pulsars are estimated. 
Subject headings: gamma rays: theory - pulsars: general 



L INTRODUCTION 

High-energy observations with the Fermi Large Area Tele- 
scope in the GeV band of four pulsars. Vela, PSR J1028- 
5819, PSR J0205H-6449, and PSR J2021H-3651 have been re- 
cently reported, and the GeV light curves of the four pul- 
sars show two main peaks separated in phase by 0.4-0.5 and 
lagged the radio peak by a phase shift of 0.08-0.2(Abdo et al. 
2009a,b,c,d). These new observed results provide important 
limitations to 7-ray emitting models of pulsars. 

Many models have been proposed to explain high-energy 
properties of 7-ray pulsars. The common feature of these 
models is that electrons and/or positrons are accelerated by 
the electric field parallel to magnetic field and produce pho- 
ton emission through different radiation processes in the open 
volume of the pulsar's magnetosphere. On the other hand, 
the main difference of the these models is that acceleration 
regions are different. For examples, photon emission occurs 
within several radii of a neutron star above a polar cap surface 
in the polar cap models (e.g., Daugherty & Harding 1982, 
1994; Sturner et al. 1995). However, in the outer gap mod- 
els (e.g., Romani & Yadigaroglu 1995; Zhang & Cheng 1997, 
1998; Zhang et al. 2004; Tang et al. 2008), both particle ac- 
celeration and photon emission take place in the outer parts of 
the pulsar's magnetosphere. 

In this paper, we focus on the interpretation of GeV light 
curves of the 7-ray pulsars which involves three dimensional 
models for 7-ray pulsars. In the polar cap models, the char- 
acteristic of the two peaks in the light curve are usually pro- 
duced when the line of sight intersects the polar cap beam, 
and a nearly aligned rotator is needed to reproduce the ob- 
served properties of the light curve (e.g., Daugherty & Hard- 
ing 1994). Based on the geometrical study of polar cap accel- 
erators, slot gap model in which particles continue to acceler- 
ate and radiate to high altitude along the last open field lines 
has been proposed (e.g., Arons 1983; Arons & Scharlemann 
1979; Muslimov & Harding 2003) and reasonably reproduces 
the Crab pulsar phase-resolved spectra (Harding et al. 2008). 
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In the outer gap and the two-pole caustic models, 7-rays are 
emitted in the acceleration gap located near the last open field 
lines, and prominent peaks in the light curve are produced by 
the caustic effects, i.e., the aberration of photon direction and 
the time-of-flight delay caused by the finite speed of light (Ro- 
mani & Yadigaroglu 1995; Cheng et al. 2000; Dyks & Rudak 
2003; Dyks et al. 2004). For a polar cap of a 7-ray pulsar, 
high-energy photons emitted along the leading field lines ex- 
tend in phase, but those along the trailing lines can be piled 
up in phase as a integrative result of the special configuration 
of the trailing lines and the caustic effects. 

In the original outer gap models, the gap starts from the 
null charge surface to the light cylinder, with the last open 
field lines as the lower boundary and a charge layer on a sur- 
face of the open field lines as the upper boundary (e.g., Cheng 
et al. 1986a,b, 2000; Romani & Yadigaroglu 1995). High en- 
ergy photons predicted by these models are produced beyond 
the null charge surface, and then only these from one pole can 
be observed. However, more recent electro-dynamical studies 
show that the gap position shifts if particles can be injected 
either at the inner boundary or at the outer one, and the gap 
can be located at the star surface if the injection rate across 
the boundary is comparable to the Goldreich- Julian value (Hi- 
rotani & Shibata 2001a,b, 2002; Hirotani, Harding & Shibata 
2003; Takata, Shibata & Hirotani 2004). Tang et al. (2008) 
proposed a revised outer-gap model to investigate the multi- 
band phase-resolved spectra for the Crab pulsar In the model, 
the gap extends significantly towards the star since a current 
carried by the pairs could be produced in the gap although 
ignoring the current injection from the inner and the outer 
boundaries. As a result, the emission from the both poles can 
contribute to the observed light curves. 

Dyks & Rudak (2003) firstly proposed the two-pole caus- 
tic model to interpret the high-energy light curves of pulsars. 
In their model, photons are emitted uniformly in the gap, and 
the gap extending from the star surface to high altitudes is 
confined to the last open field lines. Two sharp peaks with 
well-developed wings and a phase separation of 0.4-0.5 can 
be easily reproduced using the model; moreover, a bridge 
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emission and a significant off-pulse emission also appear in 
the hght curve. Therefore, the light curves of the pulsars with 
two main peaks separated by 0.4-0.5 in phase can be well 
explained using the two-pole caustic model. 

Motivated by the recent reports on the observations with 
the Fermi Large Area telescope for the four normal pulsars, 
we restudy the high-energy light curves for these pulsars ba- 
sically following the two-pole caustic model (Dyks & Rudak 
2003; Dyks et al. 2004) and find that the off-pulse emission 
level in the observed GeV hght curves compared with the re- 
sulting one with the two pole caustic model is very low, es- 
pecially, for Vela and PSR J0205H-0449. On the other hand, 
after taking the pair production process into account, Cheng 
et al. (2000) found that the distance to the null charge sur- 
face is the function of azimuthal angle (j> (also see Zhang 
& Cheng 2001; Tang et al. 2008), i.e. r^^i = rn^i(4>) and 
rnu\((l> > 0°) > i'aa\{4> = 0°), where (j) = Q° represents the mag- 
netic meridional plane. In other words, the extension of the 
outer gap along the azimuthal direction is limited by the pair 
production process. Therefore, we revise the two-pole caustic 
model as follows: high-energy photons cannot be efficiently 
produced along the field lines where the distances to the null 
charge surface are larger than /i„ 1 times of the distance 
of the hght cylinder, and we choose /in = 0.9 in this paper. 
As a result, the level of the emission of the hght curve in 
the off-pulse region can be greatly reduced, and the result is 
more consistent with the observations. We apply the model 
to the four pulsars. Vela, PSR J1028-5819, PSR J0205H-6449, 
and PSR J2021H-365 1 whose light curves obtained with Fermi 
have been recently reported, and the results of the model are 
consistent with the observations. 

2. MODEL AND RESULTS 

A retarded vacuum dipole geometry of the magnetic field 
around the pulsar is assumed in the model (see details in 
Cheng et al. 2000). Although the full three-dimensional MHD 
solutions of the pulsar magnetosphere can describe the struc- 
ture of the magnetosphere filled with charges more accurately 
(Spitkovsky 2006; Kalapotharakos & Contopoulos 2009), the 
simulations are numerical and time-consuming. The retarded 
vacuum dipole which can approximate the MHD solutions is 
analytical and easy to carry out, thus we use the retarded vac- 
uum dipole in the calculation (see also in Harding et al. 2008). 
In the original two-pole caustic model, the gap confined to the 
last open field fines is thin and extends from the polar cap to 
the hght cylinder; the high-energy photons are emitted along 
the field lines extending from the star surface to high altitudes 
and the emissivity is uniform within the gap region (Dyks & 
Rudak 2003; Dyks et al. 2004). 

The calculation basically foUows the method in Dyks et al. 
(2004). Runge-Kugga integrations are employed to receive 
the shape of the polar cap rim, and the open volume coordi- 
nates (rove, ^ovc) are established to easily emulate the particle 
distributions at the pulsar surface, where rove = 1 ± <iovc, "^ove 
is the minimum distance of a point from the polar cap region 
in units of the standard polar cap radius, /ovc is the arc length 
along the deformed ring of the fixed rove; the electron distribu- 
tion is rim-dominated at the star surface and can be expressed 
by a Gaussian function: 
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TABLE 1 

Parameters for the four pulsars 



Name 


P (ms) 


a 


c. 


fj 




Vela 


89.3 


72 


57 


0.02 


0.0425 


J2021+3651 


103.7 


70 


81 


0.02 


0.065 


J1028-5819 


91.4 


75 


65 


0.02 


0.12 


J0205+6449 


65.7 


85 


78 


0.04 


0.01 



"A phase interval 5ph is added to the resulting hght 
curve to compare with the observation for each pul- 
sar. 



ting region to the star, and the emissivity drops to zero at a 
distance p^ax = (0.75-0.95/?ic) from the rotational axis, where 
Ric is the distance of the hght cylinder (see details in Dyks et 
al. 2004). r^yg = 1 means that the particles emitting 7-rays are 
accelerated mainly along the last open field lines. However, 
we use r[Jyj, = 0.98 in the calculation as the physical outer-gap 
model in which the surface containing the last open field lines 
is usually treated as a boundary of the gap (e.g., Tang et al. 
2008). 

The high-energy photons are emitted tangentially to local 
magnetic field lines in the corotating frame, and the direction 
rj in the observer frame can be obtained from the direction 77' 
in the corotating frame with 



+ [7+ (7- ??')//?']/? 



(2) 



7(1 

where7 = (1-/3^)"^/^, /3 = v/c, v = f2 x r, fiis the angular ve- 
locity of the pulsar and r is the radial distance of the emitting 
point (Dyks & Rudak 2003). Then a phase of detection $ can 
be obtained with 
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here a describes the thickness of the gap; the emitting region 
has an upper boundary rmax. which is the distance of the emit- 



where cf)g„ is the azimuthal angle of 77. 

The resulting projection of the last open field lines on the 
space ($,Cobs) can be seen in panel (a) of Fig.l, where Cobs 
is the viewing angle. The two blank deformed ovals corre- 
spond to the two polar caps, and two caustics resulting from 
aberration and the finite speed of light can be seen as the dark 
arches along the trailing field lines in the projection map. So 
photons can be piled up along a broad range of altitude be- 
cause, for the traihng fines, the effect of different altitudes on 
the final phase can be almost compensated by the aberration 
and time of flight (Morini 1983; Dyks & Rudak 2003; Dyks 
et al. 2004), and two main peaks can be usually seen in the 
pulsar high-energy light curves. 

We now apply the model to the four normal 7-ray pul- 
sars. Vela, PSR J1028-5819, PSR J0205H-6449, and PSR 
J202 1-1-3651 whose GeV light curves obtained with the Fermi 
have been recently reported and have two main peaks sepa- 
rated by 0.4-0.5 in phase (Abdo et al. 2009a,b,c,d). Vela is a 
young pulsar with a characteristic age of 1 1 kyr and a period 
of ^ 89 ms (Large et al. 1968). The distance is ~ 287 pc from 
the VLBI parallax measurement (Dodson et al. 2003). It is the 
brightest persistent 7-ray source and was observed first during 
the SAS2 mission (Thompson et al. 1975) and followed by 
COS B (Kanbach et al. 1980), EGRET (Kanbach et al. 1994), 
AGILE (Pellizzoni et al. 2009) and Fermi (Abdo et al. 2009b). 
The Fermi telescope obtained 32,400 pulsed photons with en- 
ergies > 0.03 GeV , and the high-quality result shows that 
two main peaks, PI and P2, separated by ~ 0.432 in phase 
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Fig. 1. 

a = 12°, rmax 



(a) Projection of the last open field lines on the space Cobs)- (b) Projection of the last open field lines for which r^ui < 0.9i?ic. The parameters < 
ax = l.ORic, pmax = 0.95ffic, and P = 89.3 ms. 



1.2 
0.9 
0.6 
0.3 
0.0 
1.2 
0.9 
0.6 
0.3 
0.0 





Vela 


V — - 




— \ 1 \ 1 — 

Vela 


H ' \ ' : 



0.0 



0.2 0.4 ^ 0.6 0.8 1.0 



Fig. 2. — The light curve for Vela (dotted line) predicted by the model 
without (upper panel) and with (lower panel) the prescription that the filed 
lines with r„ui > 0.9Ric cannot contribute emission to the resulting light curve. 
The Fermi result (0.3 - 1.0 GeV) (Abdo et al. 2009b) is shown as soHd line 
for comparison. The parameters are shown in Table 1. 

appear in the light curves and a significant component in the 
bridge region (Abdo et al. 2009b). PI lags the radio pulse at 
1.4 GHz by a phase of ~ 0.13; the two peaks are asymmetric, 
i.e., P2 has a slow rise and a fast fall, whereas the fall of PI 
is slower Moreover, a third peak in the bridge component be- 
comes distinct at > 1 GeV and the peak moves to later phase 
with increasing energy (Abdo et al. 2009b). 

The resulting light curve for the pulsar with the model in 
this paper can be seen in Fig.2, and the parameters are given 
in Table 1. For Vela and the other pulsars, the viewing an- 
gle Cobs is chosen basically according to the available results 
deduced from the geometry of the pulsar wind nebula in the 
X-ray band, and then we choose the other parameters to well 
reproduce the Fermi light curves. Note that the phase for 
the LAT data corresponds to the radio peak, whereas in this 
paper it corresponds to photon emitted at the stellar center. If 
the radio pulse is emitted above the stellar surface, the phase 
of the radio peak will deviate from due to the aberration and 
retardation effects, in the resulting curve would coincident 
with the observed ones assuming the radio pulse is the core 
component. Moreover, the location of the radio emission for 
a pulsar can be 10 - 1000 km (Karastergiou & Johnston 2007; 



Watters et al. 2009) high above the stellar surface. Therefore, 
there is usually a phase difference between the resulting curve 
of the model and the observed one with the LAT for the aber- 
ration and retardation effects, and a phase interval (5ph, which 
is chosen according to the first peak in the resulting curve is 
aligned with the observed one can be obtained for each pulsar. 

Vela has a viewing angle of ^ 63° according to the torus 
fitting procedure to the Chandra observation (Ng & Romani 
2004, 2008). We use different Cobs around this value and other 
appropriate parameters to detect which one can well repro- 
duce the observed profile. The final results show that the pro- 
file of PI, P2 and the bridge emission in phase ~ 0.32-0.5 
can be well reproduced with the model and a notch follow- 
ing PI also appears in the resulting light curve. However, 
the model result underestimates the emission in the P3 region 
($ = 0.15 — 0.3) and significantly overestimates the emission 
in the off-pulse region (<1> = — 0.1,0.6—1.0). As argued in 
Abdo et al. (2009b), the emission in the P3 region is produced 
via the synchrotron emission resulting from the relatively low- 
altitude pair cascades initiated by the GeV curvature photons. 
The model in this paper assumes photons to be emitted along 
the magnetic field lines with uniform emissivity, so the model 
should not be used to tackle the P3 profile although the notch 
following PI can also be interpreted with it. 

Physically, in outer gap models when relativistic particles 
propagate outwardly along the field lines in the gap of a pul- 
sar, some other particles with the opposite charge flowing 
towards the star can emit high-energy 7-rays inwardly. In 
these models, high-energy 7-rays can be emitted by the the 
inwardly-propagating particles in the gap, and these photons 
can collide with the X-rays from the star surface although 
the inward emission is much fainter than the outward one. 
The inward emission can usually be neglected when model- 
ing the high-energy light curve of a pulsar. The inner bound- 
ary of the gap in the outer gap models is limited at the loca- 
tion where significant pairs are produced to screen the gap 
(Zhang & Cheng 1997, 1998; Cheng et al. 2000; Tang et 
al. 2008). Formerly, the inner boundary of the gap is set to 
the null charge surface (e.g., Cheng et al. 2000). However, 
more recent detailed electrodynamical studies show that the 
inner boundary can be shifted from the null surface to the 
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Fig. 3. — Radial distance to the null charge surface for the last open field 
lines, here </>p is the azimuth angle the polar cap. The parameters are the same 
as Fig.l. 

Star surface with the current through the gap taken into ac- 
count (Hirotani & Shibata 2001a; Takata, Shibata & Hirotani 
2004; Tang at al. 2008). Tang et al. (2008) used a prescrip- 
tion that jg((f))oorta,i((t>)~^^^ neglecting the out current injected 
from the inner and outer boundaries to get the location of the 
inner boundary, where jg{4>) is the current density in units of 
QB/ln carried by the pairs created in the gap and r„ui((/)) is 
the radial distance to the null charge surface on the last open 
field line for the polar cap azimuth angle (f). The pair creation 
process is sensitive to the gap geometry which is still not ac- 
curately known (Tang et al. 2008), and the characteristics of 
the current through the gap in 3-dimensional geometry are un- 
clear now. Following Tang el al. (2008), the altitude of the 
inner boundary of the gap for each field line is proportional 
to location of the null surface. So for the field lines with high 
null surface, the gap cannot be formed or effective enough 
to produce high-energy 7-rays. In fact, we found that if the 
emission from the field lines with ^nui > 0.9/?ic is excluded, 
the resulting light curve for Vela shows two peaks with zero 
flux from the pulsar magnetosphere in the off-pulse region, 
which is more consistent with the Fermi observation. The ra- 
dial distance of the null surface for each last open field Une 
is shown in Fig. 3, and the azimuth extension of the last open 
field lines emitting high-energy photons is ^ 230° for the Vela 
parameters with fin = O.9.. Therefore, different from the two- 
pole caustic model proposed by Dyks & Rudak (2003) and 
? in which the particles can emit high-energy photons along 
all the last open field Unes, the high-energy photons from the 
field lines for which rnui is larger than /i„ times of the distance 
of the light cylinder are excluded. As a result, the correspond- 
ing emission pattern of the lines is indicated in the penal (b) 
of Fig.l. 

With the prescription that high-energy 7-rays cannot be 
emitted along the lines for which rnui > 0.9R\c, a large part 
of area in the projection map is blank. In such a case, the flux 
in the light curve following the P2 is greatly reduced with 
a < 80° and just one main peak remains since the PI disap- 
pears with a < 40°; however, the profile of the resulting Ught 
curve does not change for a ^ 90° because the filed lines with 
''nui > O.9.R1C do not contribute emission to the light curve for 
this incUnation angle. The final light curve with the prescrip- 
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Fig. 4.— The light curve for PSR J2021+3651 (dotted line) predicted by 
the model without (upper panel) and with (lower panel) the prescription that 
the filed lines with ri,, > 0.9Ric cannot contribute emission to the resulting 
light curve. The Fermi result (> 0. 1 GeV) (Abdo et al. 2009c) is shown as 
solid line for comparison. Others are the same as Fig. 2. 

tion taken into account is shown in the lower panel of Fig. 2. 
The flux outside the two peaks is greatly reduced and the mod- 
eling light curve is more consistent with the observation with 
a phase interval of 0.0425 added to the resulting light curve to 
make the phase of the first peak consistent with the observed 
one. 

PSR J202 1+3651 is a young and energetic pulsar with a 
rotation period of '--^ 104 ms and was discovered by Roberts 
et al. (2002) in the radio search towards the five unidenti- 
fied ASCA X-ray sources coincident with the EGRET 7-ray 
sources. The pulsar is associated with 3EG J2021-3716 (Hart- 
man et al. 1999; Roberts et al. 2002), and the pulsed 7-rays 
were first discovered by Halpern et al. (2008) with AGILE. 
Possible pulsed X-rays from the pulsar and a pulsar wind 
nebula (PWN) associated it had been detected using Chan- 
dra (Hessels et al. 2004). The PWN has a "dragonfly" shape, 
and a viewing angle of 86 ± 1° can be estimated for the pul- 
sar from a fit to the torus structure (Van Ettan et al. 2008). 
The high-resolution 7-ray light curve of PSR J202 1+365 1 ob- 
tained with the Large Area Telescope (LAT) on Fermi has 
been reported by Abdo et al. (2009c). Two narrow peaks sep- 
arated by ~ 0.468 in phase appear in the curve, and the first 
peak lags the 2 GHz pulse by 0.162 in phase (Abdo et al. 
2009c). 

We model the GeV Ught curve of PSR J2021+365 1 with the 

parameters, P = 103.7 ms ct = 70°, Cobs = 81°, cr = 0.02 (see 
Tablel and Fig.4). For these parameters, the flux in the off- 
pulse interval has a significant amount, we choose 50 counts 
per 0.01 phase in the original observed Ught curve (Abdo et al. 
2009c) as the background emission level, i.e., the flux equal- 
ing in Fig.4 corresponds to the background emission 50 
counts per 0.01 phase. A phase interval of 0.065 is chosen 
in the figure. The GeV light curve can be well reproduced 
using the model, and, for a = 70°, Cobs = 81°, the Ught curve 
excluding the emission from the field lines with rnui > 0.9R\c 
is nearly same as that including the emission except a small 
notch around the phase 0.68. Note that the a notch also ap- 
pears in the observed light curve with the LAT, which also 
sustains the correctness of our prescription. 

The pulsar PSR J1028-5819, newly discovered in the high- 
frequency search for pulsars using the Parkes telescope and 
the Australia Telescope Compact Array on the unidentified 
EGRET source 3EG J1027-5817, is a young pulsar with a 
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Fig. 5. — Comparison of the resulting liglit curve with the Fermi observa- 
tion (> 0. 1 GeV) (Abdo et al. 2009a) for the pulsar PSR J1028-5819. Others 
are the same as Fig.2. 



PSR J0205+6449 




0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 6. — Comparison of the resulting light curve with the Fermi observa- 
tion (> 0. 1 GeV) (Abdo et al. 2009d) for the pulsar PSR J0205-H6449. Others 
are the same as Fig.2. 

period of 91.4 ms and a characteristic age of 9.21 x lO'* yr 
(Keith et al. 2008). Pulsed 7-ray signals from this pulsar have 
been recently discovered by the LAT, and two narrow peaks 
PI and P2 with a phase separation of ~ 0.46 appear in the 
light curve, i.e., PI at phase ^ 0.200 and P2 at phase ^ 0.661 
with the phase set to the 1.4 GHz radio pulse (Abdo et al. 
2009a). 

The resulting light curves predicted by the model are shown 
in Fig.5 with P = 91.4 ms, a = 75°, Cobs = 65° and a = 0.02 
with a phase interval 5^^ = 0.12 is added to the resulting light 
curve to make the two peaks consistent with the observation, 
which means the altitude of the radio emission is ~ 263 km. 
For these parameters, if all the open field Unes originating 
from the stellar surface at Tovco can contribute to the observed 
light curve, the flux at PI is high above that at P2, and a sig- 
nificant amount of emission belongs to the off -pulse interval; 
however, similar to the Vela pulsar, if excluding the emission 
from the field hues with r„i > 0.9/?ic, the flux in the off-pulse 
region is greatly reduced, and the flux at PI is smaller than that 
at P2, which is more consistent with the observation. Clearly 
from the lower panel of Fig.5, the pattern of PI, the bridge 
emission, P2 and the off-pulse interval at phase 0.6 - 1 can 
be well reproduced with the model although the model under- 
estimates the left wing of PI, i.e., in the phase interval 0.0 - 
0.16. 

The pulsar PSR J0205H-6449 has a period of 65.7 ms and 
locates in the PWN 3C58 with a distance of 3.2 kpc. Pulsed 



signals are first discovered in X-rays from the Chandra X-ray 
observatory data (Weiler & Panagia 1978; Becker et al. 1982; 
Roberts et al. 1993; Murray et al. 2002). The X-ray profile 
shows two narrow peaks separated by 0.5 in phase, and the 
first X-ray peak lags the radio peak at 2.0 GHz by ~ 0.10 
in phase (Murray et al. 2002; Livingstone et al. 2009). This 
pulsar also was recently observed with the LAT on Fermi, and 
2922 effectual pulsed photons with energies > 0.1 GeV were 
detected (Abdo et al. 2009d). The Ught curve also shows two 
peaks separated by 0.49 ± 0.02 in phase. The first peak PI 
lags the radio pulse by ~ 0.08 in phase, and the seconde peak 
P2 is asymmetric with a slow rise and a fast fall (Abdo et al. 
2009d). 

The comparison of the Fermi LAT observation and the re- 
sulting Ught curve for PSR J0205H-6449 is shown in Fig.6, and 
the parameters are P = 65.7 ms, a = 85°, Cobs = 75°, a = 0.04 
and (5ph = 0.01. For these parameters, the flux both at PI and 
in the off-pulse region (phase = 0.65 - 1 .0) is significantly re- 
duced if emission from the field lines in the emitting region 
with rnui > 0.9/?ic is excluded, which is more consistent with 
the observation. With the prescription taken into account, the 
model can broadly reproduce the GeV light curve obtained 
with the LAT for pulsar PSR J 1028-58 19, alfliough it under- 
estimates the flux in the left wing (phase = 0.0 - 0.08)(see 
Fig.6). 

3. DISCUSSION AND CONCLUSIONS 

In this paper, we investigate the high-energy 7-ray light 
curves from pulsars using a revised two-pole caustic model 
and apply it to the four pulsars, Vela, PSR J1028-5819, PSR 
J0205-I-6449, and PSR J2021-I-3651 whose light curves in the 
GeV band have been obtained with the LAT on the Fermi tele- 
scope. The calculation basically follows the two-pole caustic 
model proposed by Dyks & Rudak (2003) and Dyks et al. 
(2004). The original two-pole caustic model assumes high- 
energy particles are rim-dominated (r^^^ = 1) and emit 7-ray s 
along all the last open filed lines extending from the stellar 
surface to high altitude, i.e., the polar cap azimuth extension 
of the emitting region is 360°. However, in our revised two- 
pole caustic model, the gap accelerating particles is screened 
by the pair creation process along the azimuthal direction, and 
the last open field lines are usually treated as a boundary of the 
gap. As a result, we use r^^^ = 0.98 in this paper Moreover, 
we find that the resulting high-energy light curve for Vela is 
significantly higher than the Fermi result, and this discrep- 
ancy can be eUminate if excluding the emission from the field 
lines where rnui > 0.9/?ic. Physically, in the former outer-gap 
model (e.g., Zhang & Cheng 1997, 1998; Cheng et al. 2000), 
the inner boundary of the outer gap is the null charge sur- 
face where pairs creation process is significant; more recent 
electrodynamical studies show that if the current through the 
gap is taken into account, the inner boundary can shift to the 
star surface, and the location of the inner boundary is deter- 
mined by the current injected from the inner and outer bound- 
aries and that produced in the gap (Hirotani & Shibata 2001a; 
Takata, Shibata & Hirotani 2004; Hirotani 2006, 2008; Tang 
et al. 2008). For the field lines with relatively high null charge 
surface, the gap cannot be formed or the gap is not effective 
in producing high-energy photons. However, the details of the 
three-dimensional gap are usually uncertain for a pulsar due 
to both the sensitivity of the pair creation process in the gap 
and the unknown of the outer current. On average, it can be 
assumed the gap cannot be formed on the field lines where 
''nui > 0.9/?ic, and high-energy photons can be produced ef- 
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fectively and uniformly from the star surface to high ahitude 
along the other field lines. Note that the electric field strength 
near the inner boundary in these models is usually not big 
enough to make the local particles radiation reaction-limited 
(Hirotani 2006, 2008; Tang et al. 2008), and then the unifor- 
mity of the emissivity of the high-energy photons from the star 
surface to high altitude is violated. Therefore, our geometrical 
model for the light curves of pulsars has a discrepancy with 
the outer-gap model. 

Recently, Bai & Spitkovsky (2009a) argued that the treat- 
ment of the aberration effect with the retarded dipole mag- 
netic field in the original TPC model is not self-consistent 
because the retarded dipole magnetic field is valid in the lab 
frame rather than the instantaneously corotating frame. We 
note that the retarded dipole formula is an approximation to 
the more realistic force-free field to the first order of r/Rijc, 
while the difference between treating the retarded dipole field 
in different frames is to the second order of rjRic. Moreover, 
the calculated high-energy light curves are usually consistent 
with the observations when the retarded dipole magnetic field 
is assumed to be valid in the instantaneously corotating frame 
(see also in Tang et al. 2008; Venter et al. 2009). Therefore, in 
our paper, a retarded vacuum dipole geometry of the magnetic 
field around the pulsar is employed as in the instantaneously 
corotating frame. 

Profiles of 7-ray light curves of pulsars are sensitive to the 
structure of the magnetic field lines and the location of the 
gap, and uncertainties still exists using the vacuum magnetic 
field since the magnetic filed structure near the light cylin- 
der can be significantly influenced by plasma current (Bai (& 
Spitkovsky 2009a). Bai & Spiticovsky (2009b) modeled 7-ray 
pulsar light curves with a simulated force-free magnetic field. 
They found that double-peak pulse profiles cannot be pro- 
duced either in the conventional TPC model or in the outer- 
gap scenario using the force-free structure. Instead, an "aimu- 
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lar gap" is proposed, in which the emission zone is located at 
open field line regions that are just outside the current sheet. 
Nevertheless, the origin of the gap and particle acceleration 
mechanism jire still lacking. 

Photons with energies above several GeV can been attenu- 
ated either through 7-7 or through 7-B interaction in the mag- 
netosphere, and the profile of the light curve changes with 
different energies. As a result, the model, which assumes 
the high-energy photons are emitted uniformly along the field 
lines in the gap, cannot be used to interpret the difference of 
the light curves with different energies. On the other hand, 
the emissivity of photons with low energies about several tens 
of MeV changes greatly in the magnetosphere, and then the 
model also cannot be employed to study the MeV emission 
of pulsar. However, because the electrons can be radiation 
reaction-limited within a large range of altitudes in the gap 
(Hirotani, Harding & Shibata 2003; Dyks & Rudak 2003), 
the assumption is reasonable for the GeV photons which do 
not encountered significant attenuation. 

The four pulsars. Vela, PSR J1028-5819,PSR J0205-I-6449, 
and PSR J2021H-3651, have been observed with flie LAT on 
Fermi. The GeV light curves for the four pulsars are usually 
shown with two main peaks with phase intervals from 0.43 - 
0.5. With appropriate parameters, the GeV light curves can 
be broadly reproduced with the model. For Vela, PSR J 1028- 
5819 and PSR J2021H-3651, the emission level in flie off-pulse 
region is greatly reduced excluding the emission from the field 
lines with rnui > 0.9/?ic, which is more consistent with the ob- 
servations. 
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